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NATIONAL  ADVISORY  COMMITTEE  FOR  AERONAUTICS 


TECHNICAL  NOTE  2^36 

HEAT  DELIVERY  IN  A  COMPRESSIBLE  FLOW  AND 
APPLICATIONS  TO  HOT-WIRE  ANEMOMETRY 
By  Chan-Mou  Tchen 

SUMMARY 


In  a  ‘two-dimensional  field  a  generalized  potential  theory  applicahle 
to  nonadiabatic  and  rotational  flow  is  developed.  Three  partial  differ¬ 
ential  equations  are  first  obtained  determining  the  three  variables  which 
are:  Distribution  of  additional  temperature  -d,  velocity  perturbation 
and  an  auxiliary  function  k  characterizing  the  rotationality  of  the 
flow.  With  the  use  of  this  theory  the  action  of  heat  sources  on  the 
flow  is  studied^  and  the  heat  delivery  in  a  compressible  flow  at  subsonic 
and  supersonic  speeds  is  calculated.  The  results  show  the  effect  of 
compressibility  and  the  nonlinear  cooling.  Applications  of  the  results 
to  hot-wire  anemometry  are  discussed. 


INTRODUCTION 


In  a  two-dimensional  stationary  field,  consider  a  certain  distrib¬ 
ution  of  heat  so-urces  placed  in  a  horizontal  stream  of  main  velocity  U, 
which  can  be  subsonic  or  supersonic.  The  following  points  will  be 
studied  theoretically: 

(1)  Action  of  heat  addition  on  the  motion  of  a  con^jressible  fluid 

(2)  Effect  of  compressibility  on  the  heat  delivery 

(3)  Nonlinearity  in  the  relation  between  heat  delivery  and  temper¬ 

ature  difference 

For  the  investigations  of  these  properties,  the  following  problems  will 
enter,  and  their  analysis  will  form  the  main  part  of  the  present  paper: 

(1)  Development  of  a  generalized  potential  theory  for  the  study  of 
a  compressible  flow  which  is  rotational  and  nonadiabatic  (i.e,,  with 
heat  addition) 


2 


MCA  TN  2l^■36 


(2)  Temperature  distribution  about  a  heat  source  in  a  compressible 

flow 

(3)  Distribution  of  the  flux  of  heat  delivered  by  the  heat  sources 

(4)  Perturbation  of  the  flow  by  the  heat  sources 

This  theoretical  study  will^  on  the  one  hand^  illustrate  the  charac¬ 
teristic  phenomena  in  a  compressible  flow  with  heat  addition^  and_^  on 
the  other  hand^  it  will  find  applications  to  hot-wire  anemometry  for  the 
measurement  of  velocity  in  a  compressible  flow  at  subsonic  and  supersonic 
speeds.  In  order  to  simplify  the  calculations^  without  losing  the 
essentiality  of  the  problems^  it  will  be  assumed  that  the  heat  soirees 
are  distributed  along  a  flat  plate^  which  corresponds  to  an  infinitely 
thin  metallic  ribbon^  edgewise  to  the  flow^  so  that  no  appreciable  per¬ 
turbations  due  to  the  shape  of  the  body  w^l  be  present,  and  that  pertur¬ 
bations  are  due  only  to  the  heat  sources.^  Actually  a  cylindrical  wire 
customarily  is  used  in  hot-wire  anemometry  for  practical  reasons,  but 
from  the  present  calculations  it  will  be  seen  that  the  form  of  a  flat 
plate  is  preferable  to  the  form  of  a  blunt  body  like  a  cylinder  in  the 
hot-wire  anemometry  of  a  compressible  flow,  because  the  variations  of 
the  Mach  number  with  the  cooling  are  single-valued  for  the  former  body 
and  double -valued  for  the  latter  body.  After  the  calculations  of  the 
heat  delivery  by  heat  sources  distributed  along  a  flat  plate,  the^heat 
delivery  by  heat  sources  distributed  alorig  a  blunt  body  will  be  discussed 
also.  In  the  present  calculations  the  variations  of  the  physical  con¬ 
stants,  such  as  the  coefficients  of  specific  heat  and  the  coefficient  of 
heat  conduction,  as  well  as  the  effect  of  viscosity  and  radiation,  will 
be  left  out  of  consideration. 


This  work  was  done  at  the  National  Bureau  of  Standards  under ^ the 
sponsorship  and  with  the  financial  assistance  of  the  National  Advisoiy 
Committee  for  Aeronautics. 


^Since  the  viscosity  is  neglected,  the  terminology  "plate, "  which 
occurs  frequently  in  this  paper,  has  the  only  meaning  of  indicating 
the  location  of  the  heat  sources,  without  offering  any  perturbation  to 
the  flow,  for  example,  in  the  form  of  a  boundary  layer. 
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GENERALIZED  POTENTIAL  THEORY  FOR  ROTATIOML  MOTION  OF  A 
COMPRESSIBLE  FLOW  WITH  VARIABLE  TOTAL  ENERGY 
Fiindamental  Equations 


By  writing  the  equations  for  the  motion  of  the  fluid,  the  following 
usual  symbols  will  be  introduced: 

X,  y  coordinates  of  a  variable  point  in  two-dimensional  field 

ds  =  \/dx2  +  dy2 


vector  velocity 

components  of  vector  velocity 

pressure 

density 

absolute  temperature 

ratio  of  specific  heat  for  constant  pressure  Cp  and  for 

constant  volume  Cy  (k  =  1.4  is  taken  for  air  in  numerical 
calculations) 

gas  constant  ^Cp  -  Cy  =  Cp  — ^ 

velocity  of  sound  at  temperature  T  ^^kET  ) 
coefficient  of  heat  conduction 
local  Mach  number  (w/c) 

The  following  functions  will  also  be  introduced: 


k  -  1 

=  1  +  — - — 


a  =  yi  - 


P  = 


k 
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2h  =  pQCp/X. 


* 


h'  =  pji 


2-k 

g  =  ^h'U  =  hU 

The  above  symbols  are  frequently  used  in  the  present  calculations. 
Other  ones  will  be  introduced  wherever  needed.  The  subscript  o  indi¬ 
cates  the  state  reduced  adiabatically  to  stagnation,  and  the  subscript  » 
indicates  the  state  in  the  free  stream. 

The  total  temperatirre  at  a  variable  point  is  given  by  the  definition 


Tt 


To  +  ^(x,y) 


(1) 


where  Tq  is  a  constant  called  isentroplc  stagnation  temperature  and  ^ 

^(x,y)  is  the  additional  temperatiire  due  to  an  introduction  of  heat. 

For  the  sake  of  convenience  the  relations  between  the  different  temper- 

ature  symbols  are  shown  as  follows:  « 


Tt/To  =  1  + 
Tt/T  =  1^ 


(2) 


Also, 


T/Tq  =  Tll'l 


W' 


2CpTo 


=  t(i  -  ii-1) 


(3) 


In  the  present  problem  there  are  four  dependent  variables  w,  p, 
p,  and  tS.  For  these  four  variables  four  equations  will  be  written  by 
expressing  the  following  four  laws  for  the  motion  of  the  fluid:  ^ 

(1)  The  law  of  conservation  of  momentum 

M 
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(2)  The  law  of  conservation  of  mass 

(3)  The  law  of  conservation  of  energy 

(4)  The  law  of  state  of  a  perfect  gas 

The  conservation  of  momentum  is  expressed  by  the  equation  of 
moment\mi.  In  vectorial  form  it  can  be  written  as  follows: 


1/2  grad  w^  +  rot  w  X  w  =  -i  grad  p 

P 


In  a  two-dimensional  problem  this  equation  for  w  can  be  separated  into 
two  equations  for  u  and  v  as  follows: 


u 


Su 

Bx 


+ 


V  ^  ^ 

Sy  p  Sx 


(4a) 


i 

P  Sy 


(4b) 


The  conservation  of  mass  of  the  moving  fluid  is  expressed  by  the 
equation  of  continuity 


div  w 


w 

P  Ss 


■A 


where  ds  -  \/dx^  +  dy^  is  the  length  of  a  segment  of  the  streamline. 

^e  conservation  of  energy  is  expressed  by  the  equation  of  energy^ 
sometimes  called  the  equation  of  heat  conduction 


pw  •  grad  ^c^T  ^  j  ^  (p'W’)  +  X  div  grad  T 


This  equation  can  easily  be  transformed  into  the  following  one: 


pw  •  grad  ^CpT  +  i  w^  j  =  X  div  grad  T 


(6) 
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In  equation  (6)  T  can  te  replaced  by  t  from  definition  (l),  and 

the  following  equation  can  be  expressed:  * 

1 

p  =  Po(T/To)^-l  Si 

with 

.  (So-S)/E 

Si  =  e 


where  S  is  given  by  the  definition  of  entropy 

T  dS  =  Cp  dT  -  i  dp 


Then  equation  (6)  becomes 


div  grad  =  2hSl(T/^I)^“^  w  •  grad  t 


where  2h  =  pQ^pA* 

Finally  the  equation  of  state  of  a  perfect  gas  is 

p  =  RpT  (7) 

These  five  scalar  equations  (equations  (^a)^  (^b),  (5);  (6a)^  and 
(7))  determine  the  five  variables  in  the  problem  u,  v^  p,  p^  and  '<9. 

The  problem  is  now  to  reduce  the  nximber  of  equations  by  eliminatiug 
certain  variables.  Therefore  a  system  of  three  equations  will  be  dis¬ 
cussed  which  will  determine  the  three  principal  variables  u,  v,  and 

The  first  equation  can  be  obtained  from  the  equation  of  continuity 
(5)^  by  eliminating  p  and  p  with  the  use  of  equations  (4a)^  (^b)^ 
and  (T)  Q-s  follows; 

^ ,  w^  Sw  ,  w  St 
div  w  =  k  —  —  +  “  — 

Ss  T  OS 


4 


/ 
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or  again,  by  replacing  T  hy  ■d  from  equation  (l). 


div  w 


^  w  ^ 
bs  T  Ss 


In  the  y  coordinates  this  equation  can  he  vritten  in  its 
expanded  form  as  follows: 


.  u^\Su  uv/5v  ^  5u\  ,  L  Y^\bY  if  b'd  b'^ 


Herewith  the  first  equation  of  a  system  of  three  equations  is  obtained 
which  will  be  discussed  in  the  following  paragraphs.  The  first  equation 
will  be  called  the  equation  of  motion  of  the  fluid.  The  second  equation 
is  the  equation  of  energy  already  given  by  equation  (6),  and  the  third 
equation  can  be  derived  as  follows  from  the  condition  of  the  rotationality 
of  the  flow  which  is  given  by  the  equation  of  momentum  (4) .  Simple  cal¬ 
culations  of  the  rotational  term  in  equation  (4)  as  a  function  of  ■fl 
and  T-t  from  equation  (l)  lead  respectively  to 


rot  w  X  w 
^2 


k  -  1 


grad  log 


grad  ■S 


rot  wxw_  1  (k/1  ^-p^t 

- 72 -  =  "FTn:  \P  /P  /  -  ^  If  grad  logg  (9a) 


where  grad  loge  \^p  j  ^  J  nothing  else  than  the  increase  of  entropy 
expressed  by 


grad  logg 


grad  (S/cp) 


as  can  easily  be  proved  from  the  definition  of  the  entropy. 
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From  equation  (9a)  it  is  seen  that  an  isentropic  (S  =  Constant)  and 
adiabatic  flow  (Tt  =  Constant)  has 


rot  V  X  w 


Furthermore^  by  taking  the  rot  of  equation  (9^)^  making  use  of 
equations  (2)^  and  observing  that  rot  grad  A  =  0,  there  is  obtained 
finally 


rot  .  .  grad  X  grad  log^  rVs 


This  is  the  equation  of  rotationality  of  the  flow. 

The  two  equations  (8)  and  (lO)  determine  the  u,  v  field  as  a 
function  of  the  temperature  which,  in  its  turn,  is  determined  by  the 
equation  of  energy  (6). 

Generalized  Potential  Theory  for  a  Rotational  Flow 

A  system  of  coordinates  formed  by  the  streamlines  and  their  ortho 
gonal  trajectories  is  taken.  The  stream  function  is  used  for  the 
streamlines,  giving 


P  8y 

p  8x 


A  generalized  potential  function  9  and  an  auxiliary  function  k 
can  be  chosen  such  that 


K 

oy 


V 


K 


2 


NACA  TN  2436 


9 


or  In  vectorial  form  w  =  k  grad  9.  The  function  cp  is  called  the 
generalized"  potential  function^  because  it  comprises  the  particular 
function  for  an  "irrotational"  potential  flow  by  putting  k  =  1.  The 
generalized  potential  motion  can  now  be  determined  by  three  equations: 

The  equation  of  motion  in  the  form  of  a  differential  equation  for  cp, 
the  equation  of  energy  which  gives  the  temperature  distribution,  and  a 
third  equation  which  determines  the  auxiliary  function  k. 

It  is  worth  while  to  remark  that  the  study  of  the  distribution  of 
a  certain  physical  property  (e.g.,  temperature  distribution)  about  a 
body  of  more  complicated  form  can  be  made  possible  by  the  introduction 
of  such  a  set  of  functions  cp  and  t  in  the  generalized  potential 
theory,  because  the  use  of  such  variables  cp  and  t  as  new  coordinates 
will  transform  the  body  into  a  flat  plate. 

The  third  equation  which  determines  the  auxiliary  function  k  can 
be  constructed  as  follows.  If  the  u,  v  field  is  considered  momentarily 
as  given,  the  function  k  will  be  determined  by  the  following  relations 
derived  by  simple  vectorial  operations: 


rot  w  =  -w  X  grad  log^  k 


(13) 


From  equation  (13)  the  rotational  functions 
can  be  calculated  as  follows: 


rot  w  X  w 

c2 


and 


rot 


rot  w  X  w 

c2 


rot  w  X  w 
c2 


w2 

— 2  lose 


w  •  grad  logg  k 


V 


(14) 


.  rot  w  X  w  ^  w^ 

rot  - - -  =  -  grad  —  X  grad  log^  k  - 


w  •  grad  logp  k 


w  X  grad  log^  k 


w  X  grad 


grad  logg 


(15) 


t 
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By  comparing  equation  (15)  with  equation  (lO)^  the  equation  for  k 
is  obtained  in  the  vectorial  form  as  follows: 


w2  w  •  grad  loge  k 

grad  X  grad  loge  - 2 - ^  ^  grad  logg  k  + 

c^  ^ 

„  x  /vj_^2og^\  ,  g,,,  ^  ^  ^ad  loge  (16) 

The  velocity  potential  cp  can  he  written  in  two  parts  as  follows: 

cp=^+(i)  (it) 

where  cp  is  the  perturbed  velocity  potential,  cp  =  Ux  is  the  velocity 
potential  in  the  free  stream,  eind  $  is  the  potential  perturbation. 
Also, 


u  =  k(U  +  SiJi/Sx)  =  U  +  u' 
u'  =  K  +  (k  -  l)U  ’• 

V  =  v'  =  K  Siji/Sy 


(18) 


where  U  is  the  velocity  in  the  free  stream. 

After  substitution  of  relations  (18),  equation  (8)  becomes 


with 


(1  «  -  2  —  +  (1 

\  bxhy  \ 


c^/ 


f  = 


iL 

KT 


W 


5t 

Sb 


w  — 

/ 


(19) 


ds  =  +  dy2 

w  =  \jvL^  +  v^ 


M 


/I 


0 


vhere 
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It  may  be  useful  to  group  the  main  equations  as  follows 

(a)  Equation  of  motion  (see  equation  (19)): 


.  L 


1  -  ^  ^  -  2 


c^/hx^  Sx5y 


+  1  - 


where  f  is  given  in  equation  (19) 

(b)  Equation  of  energy: 


A(h”^T)  =  2h(T/^)^”^  w  ^ 


^2  ^2 

where  A  =  This  equation  is  the  analytical  form  of 

equation  (6a). 

(c)  Equation  of  rotationality  (see  equation  (16)): 


,  w^  ^  ^  ^  ’  gi'ad  loge  K 

grad  —  X  grad  logg  K  +  - - - w  X  grad  logg  k  + 


w  X  grad 


w  •  grad  loge  k 


grad  X  grad  logg 


If  s  and  n  represent  the  unit  vectors  along  a  streajnline  and 
along  a  normal,  respectively,  the  following  equations  can  be  written: 


w  •  grad  loge  k  =  w 


S  loge 


w  X  grad  logp  K  =  w 


S  loge  K 


j  >.2  V.  on  ^  loge  ^  loge  k 

grad  X  grad  logg  k  =  - -2H - iii - — 

os  on  on  3s 
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Using  these  formulas  and  assuming 


S  log^  K  h  log^  K 
- 2 — 

hs  on 


A 


(20) 


there  is  obtained  from  equation  (l6) 


K  « 


(21) 


In  an  adiabatic  flow  (i3  =  O),  the  expression  (2l)  reduces  to  the  well- 
known  value  K  =  1.  Assumption  (20)  expresses  the  condition  that  the 
gradient  of  the  additional  temperature  (or  additional  energy)  must  be 
much  smaller  along  the  streamlines  than  along  their  orthogonal  trajec¬ 
tories.  This  condition  exists  in  problems  of  the  introduction  of  a  heat 
boundaiy  (problems  of  combustion)^  of  heat  sources^  or  of  a  hot  body  in 
a  stream  moving  at  a  high  speed.  In  all  those  problems  the  motion  of  the 
fluid  is  extraneous  to  the  heat  introduction*  Othervise  condition  (20) 
will  not  be  valid,  for  exanple,  in  the  case  of  an  explosion,  where  the 
predominant  motion  of  the  fluid  is  created  by  the  heat  emission  Itself. 


In  the  present  paper  the  problem  of  the  rotational  and  nonadiabatic 
flow  will  be  treated  under  the  assumption  (20)  and  hence  with  the  value 
of  K  given  by  expression  (21).  Siibstltuting  expression  (21)  into 
equation  (19),  the  equation  of  motion  becomes 


with 


/  u^\S££  _  2  —  +  (l 

\  c^/hx^  c^  hxby  \ 


c2jay2 


f  =  ^(l  +  kM^)T"3/2^^ 


(22) 


Small  Perturbations 

As  a  first  approximation,  equations  (6a)  and  (22)  can  be  transformed 
into  linear  ones  by  introducing  the  assumption  of  small  perturbations; 
that  is,  it  is  assumed  that  the  perturbations  u'/u,  v'/U,  $/cp,  and 

i3/Tq  flnri  their  differentials  are  small  compared  with  unity,  such  that 

terms  of  second  order  are  negligible. 
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Under  this  ass\in5)tion,  the  equation  of  energy  (6a)  and  the  equation 
of  motion  (22)  become,  respectively, 

(a)  Equation  of  energy: 


where 


^  ^  =  P  ^ 

Sy^  ^  hx 

2-k 

'4-1 

g  =  Hoo  hU 


(23) 


(b)  Equation  of  motion: 


where 


2  ^  ^  ^ 
Sx2  Sy2 

a  =  (l  - 


(24) 


DISTRIBUTION  OF  TEMPERATURE  AND  HEAT  FLUX  FOR 
A  HEATED  BODY  IN  A  COMPRESSIBLE  FLOW 
Temperatiire  Distribution 

The  convection  of  heat  is  governed  by  equation  (23).  It  has  the 
following  solution: 


^  =  Ae®^Q(gr) 


(25) 


Ik 


NACA  TW  2436 


A 

where  r  =  and  Kq  is  a  Bessel  function.  For  larger  values 

of  gr,  Kq  can  be  expanded  into  a  series  which  will  be  given  later  in 

this  section.  The  constant  A  can  be  determined  if  the  heat  delivery  Q 
is  known.  Consider  first  the  case  of  a  line  source.  The  formiola  of 
temperature  distribution  about  a  finite  body  can  be  constructed  by 
considering  a  continuous  distribution  of  heat  sovirces  on  the  body,  and 
will  be  studied  in  the  following  section. 

The  heat  delivery  Q  should  be  equal  to  the  total  heat  flux  around 
a  closed  circuit  formed,  for  example,  by  a  circle  of  radius  r  enclosing 
the  line  source.  The  following  formoila  can  be  written: 


Q  = 


+  pCpU'S  cos 


(26) 


where  ds  is  a  segment  of  arc  on  the  circle  and  r  and  tj  are  respec¬ 
tively  the  radius  and  the  angle  in  the  polar  coordinates. 

The  question  is  in  what  sector  of  the  circuit  the  heat  flow  will  be 
effective.  At  large  velocities,  the  major  part  of  the  heat  flow  occtus 
through  a  very  narrow  sector  so  that  the  magnitude  Q  will  be  at  first 
approximation  independent  of  the  limits  of  Integration,  provided  these 
are  taken  of  sufficiently  large  range.  Therefore  as  a  first  approxi¬ 
mation  there  can  be  written  in  place  of  formula  (26),  for  all  kinds  of 
velocities  (subsonic  and  supersonic). 


jt 


Q  =  2r\  j  +  2gPoo‘^-i3  cos  t\ 


(26a) 


Putting 


In(z)  =  r  /  dTje 


z  cos  n 


ii(^)  =  7  / 

'-'0 


Z  cos  T) 


cos  T) 


(27) 


i 
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formula  (26a)  can  be  written  as  foUows: 

Q  =  2AAgr  -Ki(gr)  F  drie®^  ^  + 

L  ^0 


(24„-1  -  l)Ko(gr)  r  dne^  1 


'0 


cos  T] 


=  2nXAgr  J 


-Ki(gr)lo(gr)  +  KqIi 

-  2(1  -  m«,-i)koIi  . 

By  the  use  of  the  following  expansions. 


Ko(z) 


fn  e 


-z 


\/s(^  • 


“  0  *  s) 


ii(2)  -  ^  .  2.] 

\J^\  &z) 


formula  (28)  becomes 


Q  =  2jtXAp„"^ 

HerefrcM  A  Is  obtained,  and  Horn  fonmla  (25)  there  Is  obtained 


■  sir 


23tX 


(28) 


(29) 
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Flux  Distribution  and  Heat  Delivery 

The  distribution  of  temperature  about  a  line  source  is  given  by 

^  (oQ\  This  formula  is  valid  for  subsonic  and  supersonic  flows, 

formula  {29).  This  lormuia  is  nfp'idf  represent  the  heat 

consider  nov  a  flat  plate  of  length  la  «J''^^T^nUeratnre 

flux  from  an  elementary  segment  d|  of  the  piaxe.  h> 

contribution  by  this  segment  at  a  variable  point  x,y  i 


d^  =  ^  d|q(l)eS^^-^)Ko(g\/z2  +  y2 


where  z  =  x  -  I-  By  integrating. 


g(x-l) 


kJr\Iz^  + 


Formula  (30)  gives  the  dlstrlhutlon  of 

fiS,  in"?Jr=aI^^a?i^  which  follow  fo™^a  (30)  by  the  integral 
equation  (3l)  for  the  surface  condition  on  the  plate. 

On  the  surface  of  the  plate  there  can  be  written  for  0  ^  x  ^  4a 


h  tt\ 


|g(x  -  6^ 


-g( l-x) 


[g(l  -  x]| 


or^  by  using 


the  asymptotic  expansion  of  K^{z)  for  large  values  of  z, 


d|q(  I) 


^x  -  |)  2«\ 


ft\  -2s(^-x)  r~ 

+  \2g(|  -  X 


(31a) 
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The  second  integral  on  the  right-hand  side  decreases  rapidly  with 
and  therefore  as  a  first  approximation  it  can  he  written  as  follows: 


-2g(^x)  I - S - 

d|q(|)e  \ hnZ - T 

V2g(l  -  x) 


«!■ 


2g^l  1 


For  the  sake  of  simplification  of  writing,  introduce 


-  i7  eS: 


where 


g*  =  ^2g/it 


(31b) 


The  relation  (31a)  between  the  heat  flux  and  the  wall  ten^jerature  then 
becomes 


:c=?  '3i(x) 


=  qi(x)  +  g^ 


qi(^) 


X  -  f 


First  the  integral  equation  (32)  will  be  solved  for  the  general  case  of 
a  variable  wall  temperature  Later  the  particular  case  will  be 

derived  of  a  constant  i3^  which  is  true,  if  the  plate  is  supposed  to 
be  a  very  good  conductor*. 
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The  integral  equation  will  be  reduced  to  a  linear  differential 
equation  of  first  order  by  a  differentiation  with  respect  to  x.  Firsts 

put  X  -  I  =  in  order  to  prevent  the  singularity  which  would  occur 
when  I  =  X.  The  differential  with  respect  to  x  of  the  Integral  in 
the  right-hand  member  of  equation  (32)  is 


qpd) 


q3_(x  - 


d^qi(x  -  ^2) 


qi(0)  px  q^d) 

-  +  /  dl  = 

\/x  Jq  - 


where  qj_(l)  = 

By  applying  this  rule  to  the  differentiation  of  equation  (32),  there 
is  obtained 


•  f  ^  ^ 


X  qpd) 

df  — -  =  e  •3 


This  formula  is  fulfilled  by  the  flux  at  the  distance  x.  The  flux 
at  the  distance  Xg  will  fulfill 


qi(X2) 


qi(o) 


^  in  (5) 


fXo  -  I 


=  Mw(^2) 


0 
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Multiply  every  term  of  equation  (34)  by  dx2g*(x  - 
integrate  with  respect  to  X2  between  the  limits  0  and  x,  and  sub¬ 
tract  from  equation  (33);  it  follows: 


qi(x) 


+  g* 


qi(o) 

(X 


(35) 

The  first  integral  on  the  left-hand  side  of  this  equation  has  the 
value  n.  In  the  double  integral  the  order  of  the  integration  can  be 
inverted  as  follows: 


The  last  integral  with  respect  to  X2  has  the  value  jr.  Hence  the 
value  of  the  double  integral  becomes 
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with 


(38) 
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In  the  particular  case  of  a  plate  vhich  is  a  good  conductor,  the 
vail  temperature  can  he  assumed  constant;  then  formula  (37)  becomes 
simply 


qi(x)  =  q^(0)e^®^^l  -  erf  \/2gx  j  (37a) 

or^  by  retiirning  to 

q(x)  =  q(0)e^^^^l  -  erf  ^2gxj 

This  relation  is  plotted  in  figure  1,  where  ^gx  is  denoted  by 
which  can  be  expressed  as 


u 


= 


2-k 

2(k-l) 


where  P  is  the  Peclet  number  as  will  be  defined  later. 

The  constant  qp(0)  in  formula  (37^)  can  be  determined  by  putting 
X  =  0  in  equation  (32)  which  gives: 

=  2^  =  SA  (37b) 

Now  pass  over  to  the  calculation  of  the  heat  delivery.  By  defi3altion^ 


Ql/qi(0)  =  ^  dxe^®''  erf  ^ 
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The  first  integral  of  the  right-hand  member  is  evident.  By  changing 

P 

the  variables  t]  =  2gx^  the  second  integral  of  the  right-hand  member 
can  be  written  as  follows: 


Hence^ 


Q 


qi(o)  2g 


ySga  -  l]  +  -  erf  ^/Bga  ) 


-ji  (39) 


where^  for  ga  »  1^  Oj^  is 


cr-,  =  —  \/8ga 


(39a) 


Define  a  Nusselt  number  N  such  that 


Q  =  laA  33v/4a 


vhere  A  is  the  area  per  unit  length  8a.  Hence, 


N  = 


Q 


2X^3. 


w 


m 
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From  expression  (3Tti)^ 


q(0)=2W^^  (41) 

H-oo 


Hence^ 

^  ^  =  f  \/^/^^co  (^2) 


The  above  calculations  give  a  picture  of  the  distribution  of  the 
heat  flux.  However,  in  the  order  of  approximation  adopted,  it  seems 
that  the  pattern  of  distribution  of  the  heat  flux  does  not  influence 
very  much  the  heat  delivery.  This  is  readily  seen,  since  the  expression 
for  the  heat  delivery  (formula  (42))  can  also  be  obtained  by  the 
following  more  approximate  method  which  does  not  take  into  account  the 
details  of  the  flux  distribution.  To  this  end,  if  in  equation  (31) 
a  mean  heat  flux  q  is  introduced  and  x  is  put  equal  to  4a, 
the  following  equation  can  be  written: 


d?ee(^^-^^K^[g(4a  - 


due'^KQ(u) 


(42a) 


This  integral  is  known  and  has  the  value  \jQnga  -  1.  Then, 


q  =  2k'8^ngu„,“^  ^\/8jtga  -  l^"^ 
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and  N  which  Is  defined  by 


4aq 


■becomes 


N  =  +  1)  (^3) 

By  comparing  formulas  (42)  and.  (43)^  it  is  seen  that  formula  (43)  gives 
a  very  good  approximation  to  formula  (,A2)  which  was  obtained  from  the 
considerations  of  the  exact  distribution  of  heat  flux. 

A  generalized  potential  theory  has  been  developed  and  three  differ¬ 
ential  equations  have  been  obtained  which  determine  the  three  variables: 
The  distribution  of  temperattire  difference  S,  the  potential  pertur¬ 
bation  and  the  auxiliary  fimctlon  of  rotatlonality  k.  As  a  first 
approximation^  the  equation  has  been  linearized  by  neglecting  terms  of 
the  second  order  in  pert\urbations  u'/U^  '^/'^o  (see  section 

Small  Pertiurbations ) .  Then  the  equations  can  be  solved  for  ( section 
Temperature  Distribution)^  the  heat  flux  can  be  calculated  (this  section)^ 
and  hence  the  heat  delivery  can  be  obtained  as  in  formula  (43).  How¬ 
ever,  the  results  show  only  a  linear  relation  between  the  heat  delivery 
and  the  temperature  difference.  In  order  to  show  the  nonlinear  char¬ 
acter,  a  second  approximation  must  be  pursued  by  keeping  the  second- 
order  terms  of  perturbations  in  the  calculations  of  heat  delivery.  To 
this  end,  the  distribution  of  temperatirre  difference  S  obtained  in 
the  first  approximation  will  be  used  to  calculate  first  the  potential 
perturbation  $  (two  following  sections),  and  herefrom  the  velocity  per¬ 
turbation  u'/U.  An  "effective  velocity"  which  is  the  average  perturbed 
velocity  of  the  fluid  along  the  plate  can  be  Introduced  in  the  expression 
of  the  heat  delivery  which  will  then  be  of  nonlinear  nature  ( see  first 
section  under  Nonlinear  Cooling  of  a  Heated  Body  at  Supersonic  Speed) . 

It  is  to  be  remarked  that  the  process  mentioned  above  will  give  a 
second  approximation  with  sufficient  accuracy  without  entering  into  the 
Integration  of  the  nonlinear  partial  differential  equation  in  the  form 
of  equation  (6a). 
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ACTION  OF  HEAT  ADDITION  ON  MOTION 
OF  A  COMPRESSIBLE  FLOW 
Flow  Perturbations  by  Heat  Sources 
at  Subsonic  Speed 

In  this  section  the  flow  pert\irbation  at  a  subsonic  speed  by  heat 
sources  distributed  along  a  flat  plate  will  be  investigated  by  using 
equation  ( 24) .  By  means  of  the  transformation 

— 

=  X 

y^  =  ay  ►  (L4; 

=  v$ 

equation  (24)  can  be  brought  to  a  Poisson's  equation  of  the  form 


- -  +  - -  =  fVa"2 


^Xi^  Sy  2 


For  the  sake  of  mathematical  convenience  this  equation  for  the 
motion  of  the  fluid  about  a  plate  can  again  be  transformed  by  conformal 
representation  into  an  equation  for  the  motion  of  the  fltiid  about  a 
cylinder  of  radius  a,  using  the  formula  of  conformal  transformation 

z  =  Z  +  - 

where  z  =  x^  +  iy^^  represents  the  plane  of  the  flow  about  the  plate, 

and  Z  =  X  +  lY  the  plane  of  the  flow  about  the  cylinder.  By  trans¬ 
forming  equation  (45)  into  an  equation  for  a  cylinder,  there  is  obtained 
the  following  equation: 


5 = 
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This  equation  can  be  made  dimensionless  by  dividing  X  and  by 
a,  and  by  2Ua.  When  R  and  r\  are  taken  as  the  polar  coordinates 

of  the  Z -plane,  and 


dz 

dZ 


2 


=  1 


2R’’^  cos  2t]  +  R 


-4 


the  right-hand  member  of  equation  (46)  becomes 

fi  =  i  va-2(l  +  -  2R-2  cos  2ti  + 


where  3i3/Sx  should  he  obtained  hy  differentiating  expression  (30)  for 
the  temperature  distribution.  However,  the  latter  egression's  being 
of  integral  form  will  make  the  Poisson's  equation  (46)  rather  cumbersome 
to  handle.  As  a  convenient  but  not  as  a  necessaiy  step,  the  perturbation 
function  (I>,-  is  considered  as  a  superposition  of  elementary  perturba¬ 
tions  corresponding  to  elementary  heat  sources.  First  calculate  the 
elementary  perturbation  produced  by  an  elementary  heat  source  from 
equation  (46).  An  integration  according  to  the  proper  distribution  of 
heat  sources  will  give  us  the  resultant  perturbation. 

Following  the  above  procedure,  consider  first  a  heat  source  of 
intensity  Q,  placed  at  the  origin  of  the  coordinate  (center  of  the 
plate);  then  6d/6x  can  be  calculated  from  formula  (29),  and  expres¬ 
sion  (4t)  can  be  written  as  follows: 


f^  B^v 


k=  .  xi/ri) 


\/gari 


dz 

dZ 


(47a) 


where 


a 


-2 


(l  +  kM^2j, 


^w  Q 
To 


ga 


It  is  worth  while  to  transform  this  formula  into  the  polar  coor¬ 
dinates  R  and  T)  of  the  Z-plane.  To  this  end,  note  first  that  the 
exponential  term  in  the  function  f^  can  be  reduced  to  the  form 

-constant  X  ga(R-l)^ 

fi  e 
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•\ 


n 


Since  large  values  of  ga  must  be  dealt  with,  the  function  f^  keeps 
a  significant  value  only  for  ^  =  R  1  «  1.  Therefore  the  approxi¬ 
mation  ^  «  1  can  be  introduced  in  the  transformation  process.  Under 
such  an  approximation,  the  variables  x^  and  in  terms  of  ^  and 

T)  become 

x^  =  2  cos  TJ 
y^  =  2^  sin  T] 

with 

^i  = 


It  follows  that,  at  a  variable  point  downstream  of  the  source  for 
-Jt/2  «  Ti  «  +  Jt/2, 


fi 


2. 

B.v  siagt]  2 

\^gi~coFli\cos  rij  ^  cos  tj  ^ 


(i^Tb) 


Now  the  integration  of  Poisson's  equation  (46)  for  the  Z-plane  will 
be  considered.  It  is  to  be  remarked  that  the  effect  of  the  presence  of 
a  right-hand  member  in  equation  (46)  is  equivalent  to  that  of  a  certain 
continuous  distribution  of  sources  throughout  the  region  of  flow,  the 
intensity  dl  of  the  source  distribution  in  an  element  of  area  Z  dZ  dX 
being 


dl  =  ^  Z  dZ  dX 


where  Z  and  x  are  the  polar  coordinates.  The  potential  function 
induced  by  this  distribution  of  sources  can  be  calculated  by  the  method 
of  image  (references  1  to  3) • 

Consider  first  a  imit  source  placed  at  a  point  Zq  =  Ze^^  in  the 
flow  field  outside  the  cylinder.  The  velocity  potential  induced  by  this 

sotirce  at  a  variable  point  Z  =  Re  ^  is  equivalent  to  the  velocity 
potential  of  the  unit  source  placed  at  Zq,  its  imit  image  placed  at  the 

point  Zp  =  Z"^e^^  reciprocal  to  Zq,  and  a  unit  source  of  intensity  -1 
placed  at  the  center  of  the  cylinder.  When  the  potential  function  and 
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stream  function  of  the  perturbation  induced  by  the  unit  source  are 
denoted  by  and  respectively^  there  can  be  written 

0*  +  it^  =  loge  -  Zq^  +  loge  ^Z  -  -  loge  Z 

By  taking  the  real  part  of  the  above  complex  equation^  there  is  obtained 

$*  =  I  loge  -  2ZR  cos  ( X  -  Ti)  +  r2  j  + 

|[i  -  2(ZR)-^  cos  (X  -  Tj)  +  (ZR)-2j 

Evidently  this  expression  can  be  written  in  the  two  following  identical 
forms : 


$*  =  I  loge  2^  +  1  lose  [l  -  2li  cos  (x  -  tj)  +  + 

i  loge  [l  "  2I2  '^o®  (X  -  n)  +  ^2^ 

=  I  loge  +  I  loge  [1  -  2I1'  cos  (X  -  tj)  +  + 

I  loge  [1  -  2I2  cos  (X  -  n)  +  12^  (^) 


where^  for  the  sake  of  simplification^ 

ll  =  R/Z 
ll'  =  Z/R 
I2  =  (ZR)-l 
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The  "total  induced  velocity  potential  can  be  obtained  by  integrating 
for  the  whole  region  outside  the  cylinder  as  follows: 


®i 


dX2fi(Z,X)®*(  Z,X) 


(49) 


Thus  the  solution  of  Poisson  eqTiation  (46)  is  reduced  to  the  solution 
of  the  double  integral  (49). 

Substituting  expression  (48)  into  fomnila  (49),  the  double  integral 
can  be  written  as  follows: 


®i  =  f  l°ge  +  I  lose  [l  -  2li'  cos  (x  -  t])  + 

•*-  I  loge  jl  -  2I2  cos  (X  -  Tj)  +  |2^ 

^  n“  o2rt  r 

hj  J  lose  +  I  loge  (1  -  2li  cos  (x  -  n) 

^1^]  i  [1  -  2I2  cos  (X  -  T))  + 


Only  the  velocity  perturbation  at  the  surface  of  the  cylinder  (R  =  l) 
is  of  interest;  this  is 


with  I  =  l/z. 

It  can  he  proved  that  for  a  flat  plate  the  perturbation  upstream 
of  the  heat  source  decreases  very  rapidly  with  the  distance^  in  the 
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proportion  and  hence  it  is  negligible  with  respect  to  the 

perturbation  downstream.  This  means  that  the  variations  of 

have  to  be  restricted  to  the  half  plane  corresponding  to  -it/2  and 
jt/2.  Therefore  the  limits  of  integration  in  formula  (50a)  will  be 
-j(/2  and  n/2.  In  order  to  carry  out  the  integrations  of  the  right- 
hand  member  of  formula  (50a),  the  following  sing)lifications  can  be 
introduced: 

(1)  The  variable  Z  can  be  changed  into  1  +  ^  with  t,  «  1 

(2)  Since  logg  Z  a*  the  first  double  integral  of  the  right-hand 
member  of  formula  (50a)  can  be  neglected  in  conparlson  with  the  second 
one 

(3)  The  order  of  integrations  in  the  second  double  integral  of  the 
right-hand  member  can  be  interchanged^  by  first  integrating  f^  with 
respect  to 

d5fi(Z,X)  =  B.V  (ga)"^tgX 

2V2 

Hence  fomnila  (509-)  is  reduced  to 


sin^X  1  1  -  cos  (X  -  Ti) 

dX - loge  - - - - 

COS  X  sin  X  1  -  cos  (X  +  ti) 
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Now  the  integral  on  the  right-hand  side  has  the  following  limit  values 

(1)  For  T)  =  0  it  is  zero 

(2)  For  11=^  it  is 


I 


«/2 


dX 


log 


sin  X 


sin  X  1  +  sin  X 


(See  reference  4.)  For  intermediary  values  of  T),  the  approximate 


form  sin  t)  is  suggested.  Hence^ 


$ . 

10 


-B.v  ^  (ga)-2 


6\[2 


cos  T] 


(51) 


The  distribution  of  the  potential  perturbation  on  the  compres¬ 

sible  field  along  a  plate  of  length  4a  can  be  obtained  by  multiplying 
formiila  (51)  by  2Ua,  by  using  relations  (44),  and  by  replacing 
cos  T)'  by  x/2a;  then. 


®o  =  -Bl  0  (ga)-2ua  f:[l  -  (x/2a)2]3/2 
or,  by  substituting  for  Bj, 

Where  Bg  represents  for  the  sake  of  simplification 

Bg  =  i  a-2(l  +  kM„2)  ^ 


(52) 


(53) 


Formula  (52)  gives  the  potential  perturbation  by  a  single  heat  source  at 
a  point  X  from  the  source. 


32 


NACA  TN  2^36 


Now  consider  a  continuous  distribution  of  heat  sources  such  that 
the  intensity  in  the  region  between  ^  and  &  +  d^  from  the  leading 
edge  of  the  plate  is  q(l)  d|.  Thus  the  elementary  perturbations  by 
q(|)  d|  along  the  downstream  part  of  the  segment  of  plate  which  has  a 
length  4a  -  I  can  be  obtained  from  formula  (52)  by  replacing  2a  by 
4a  -  I  and  x  by  x’  -  where  x*  is  the  distance  from  the  leading 
edge.  For  x'  -  1^0, 


d$o( I) 


=  -f  B2Uh^ 

o 


q(l)  d|  -  I 

2X^3^  X*  -  I 


The  total  downstream  perturbation  corresponding  to  a  continuous 
distribution  of  heat  sources  along  the  whole  plate  is  obtained  by  an 
integration  of  the  above  expression  with  respect  to  |  between  the 
limits  0  and  x*  as  follows: 


jr  ^^cx) 

f  BpU  — 
4  2  2e: 


2g 


cl?q(l) 


4a  -  I 

X'  -  I 


X'  - 

4a  -  I 


In  order  to  calculate  the  Integral  in  the  right-hand  member  of  the 
above  expression,  the  following  simplifications  are  introduced: 

(l)  The  asymptotic  form  of  the  formula  (37a)  for  flux  distribution 
will  be  used 


where 


q(l)  ««  q(o) 

pngl. 


q(0)  = 


f—  — 

\2g  2X\I 


(2)  A  change  of  variables 
into  the  following  type: 


^  =  cos  u  transforms  expression  (5^) 
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where  F(n)  is  a  trigonometrical  function  which  appears  in  the  course  of 
the  transformation  and  v  is  a  parameter  which  turns  out  to  be 
V  =  x'/i)-a.  The  detailed  form  of  F(ti)  will  not  he  written  out  here. 

As  the  fractional  term  has  predominant  values  in  the  proximity  of 
cos  u  {»  V,  there  can  he  written  hy  approximation 

sin  \  „/  \  sin 

/  —  F(t))  r#  F(v)  /  dTi  '  — 

^cos-ly  Vv  -  cos  u  ^cos-ly  Vv  -  cos  u 


The  integral  on  the  right-hand  side  can  he  immediately  calculated^  and 
there  is  obtained 


Ua  -1/2,,  w,  >3/2 
V  (1  -  v)(l  +  v)-*/ 

VSitga 


Hia  gradient  at  the  center  of  the  plate  will  he  taken  as  the  effective 
gradient  e  which  is 


6 


1  ^ 
U  dx 


_  3\/3 


\ptB2 


=  M 

16 


\/Ka“^(l  +  kM  2) 

m 


(55) 


It  is  seen  that  in  a  subsonic  field  the  fluid  increases  its  velocity  by 
passing  over  the  heat  sources. 


Flow  Perturbations  by  Heat  Sources  at  Supersonic  Speed 

In  this  section  the  flow  perturbations  at  a  supersonic  speed  by 
heat  sources  distributed  along  a  flat  plate  will  be  investigated  by  using 
again  equation  (2h).  The  equation  will  be  first  solved  for  a  single 
heat  source  placed  at  the  origin.  The  perturbations  of  the  flow  about 
a  flat  plate  will  be  derived  by  a  superposition  of  heat  sources  as  done 
in  the  preceding  section.  Equation  (24)  can  be  made  dimensionless  by 
dividing  the  coordinates  by  a,  and  the  potential  perturbation  by  2Ua, 
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where  ij-a  is  the  length  of  the  plate  and  U  is  the  free-stream  velocity. 
Thus, 


where 


2^  ^  ^ 


1  . 

f .  i  (1 .  — 

r  =  (m2  -  l)-l/2  =  tgTi, 


(56) 


t]q  =  Mach,  angle 


By  means  of  a  change  of  variables 


Xf  =  ax 


yn  =  y 


(57) 


equation  (56)  can  he  transformed  into  the  following  form: 


=  f 

Sxi2 


(56a) 


Equation  (56a)  has  its  characteristic  lines  given  by  yj_  ±  Xj.  =  C;  these 

are  two  families  of  straight  lines  making  an  angle  of  45°  with  the 
x^-axis. 

The  problem  is  to  find  a  solution  of  equation  (56a)  such  that  it 
satisfies  the  following  two  boundary  conditions: 

(a)  The  term  S^/Syp  vanishes  at  yi  =  0 

(b)  The  perturbation  $  vanishes  along  the  front  shock  waves,  which 
are  formed  by  the  characteristic  lines  originating  at  the  leading  edge. 
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In  particular,  the  solution  at  various  distanc 
of  special  interest. 


es  on  the  x^-axis  is 


ss-  si'-.ri.'i'S-srs-;,:: 

a  line  x^p  parallel  to  Oti'  (see  fie  P)  Tn  +>,^^ 

.  ^  I  Vi=ee  iig.  In  the  plane  x.y.  apply 

the  Gaussian  theorem  in  the  form 


hV  _  ^ 

Sxi  “  Syi 


dxi  ^  ^  = 


(U  dxi  +  V  dyi) 


where  the  double  integration  is  extended  over  the  area  Ox.p,  and  the 

If  in  equation  there  arnn^Lte^  aafT^ 

follows  from  equations  (56a)  and  (58)  that 


^  dx.  +H 


=-JJ  (59) 

Ox^p 

It  IS  evident  that  .  dy,  on  the  line  Op,  and  dx^  =  -dy,  on  the 

line  Xj^p.  Consequently^  ^ 


m 


3*1  "  ayi  "■‘l)  -  S  -  »o 


"0 


r 

L±  ^  ^Ji 


—  dxi  t  —  dj-ij  .  -Oj,  1  fe 


0  VS>^i  3yi  V  I  "Syi 
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and  it  follows  from  equation  (59)  that 


Xi  ,  6$ 

-  2$p  +  'I'o  +  j 


Now,  from  the  boundary  conditions, 


gl-=_  dXidYif 


Ox^p 


=  0 


yi=o 


•p  =  ^o  =  0 


Hence, 


®^i  ~  ~  II  ‘^^i  ‘^yi^ 

OXiP 

The  above  formula  can  be  transformed  from  the  Xj^,yj^  plane  into 
x,y-plane  according  to  relations  (57).  Omitting  the  details,  the  result 
is 


2  , 

tx--5  / 


px/ 2  0^1 

dxi  /  dyif(Xiyi) 
^0 


(61) 


where  l/a  is  the  ratio  between  the  two  areas  of  integration  in  the 
respec'tive  planes 


1  ^  6(x,y) 

6(xi,yi) 
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For  a  heat  source  placed  at  the  origin,  f  is  given  by  expres¬ 
sions  (56)  and  (29),  and  there  results 


„  „  ,  1  -gar(l-x/r) 

^  ^  (1  -  x/r) 


with 


k\J^  ^  ^oc,^)^"oo  2X-d„ 


o - w 


(62) 


It  is  to  be  noted  that  an  integration  of  f  with  respect  to 
means  an  integration  with  respect  to 

V  =  -  x) 

while  keeping  x  constant.  After  some  simplifications. 


dyif(xi,yi) 


B- 


(ga) 


■2  2^ 


,\jga{\j2-l)x2^ 


dw^e”^ 


0 


By  substituting  this  into  formula  (61),  it  follows  that 


-  -o  lr„„^-2  r^'  ■  1 


X  a 


0  ^  *0 


/2  dx.  p\jga{\f2-l)x^ 


2  “V^ 

dw  e 


and,  introducing  a  new  variable  u  such  that 


u  =  \lga{\f2  -  l)x- 
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Hence  formula  (64)  becomes 


(B  ^1  /  \  -2  1  t — 

®x  —  “7^  (S3.)  logg  J3^  -  Constant 

a  \/2  ^  ^ 

As  the  constant  does  not  play  a  material  part  in  the  velocity  pertuV' 
bation,  it  can  be  disregarded,  and  there  is  obtained 


^  (ea)-  log,^ 


and^  by  substituting  for  Bj*  from  expression  (62), 


\  =  -i  (1  +  l«„2)lajga)-l  ^  ^  log^  ^ 


Now  consider  a  continuous  distribution  of  heat  sources  over  the 
length  4a  of  the  plate,  such  that  the  intensity  of  the  source  on  the 
se^ent  between  I  and  |  +  d|  is  q(|)d|,  where  q  has  the  distrib- 
u  ion  law  given  by  formula  (37a).  The  elementary  potential  pertur¬ 
bation  d3>x  due  to  q(|)  d|  at  a  point  situated  at  a  distance  x  from 

the  leading  edge  on  the  x-axis  can  be  easily  obtained  from  formula  (65) 

y  changing  in  formula  (65)  Q  into  q(5 )  d|  and  x  into  x  -  f . 

Then,  ^ 


\(2  -  1 


(x  -  I) 


.  perturbation  upstream  of  the  source  is  again  negligible  as 

in  the  subsonic  case,  the  perturbation  corresponding  to  a  contin¬ 
uous  distribution  of  the  heat  sources  along  the  total  length  of  the 
plate  can  be  obtained  by  integrating  the  above  expression  between  the 
limits  0  and  x  of  the  variable  i.  Thus, 


NACA  TN  2^36 


=  -"T  —  /  d| 

lOCT  '*'0  *^0 


—  /^  d{  SLil  log  ga  -  (x  -  I) 

To  I  L  =  J 


For  the  sake  of  simplification  in  the  following  calculations,  the 
integral  in  formula  (66)  will  he  represented  hy 


S2=  sa 


^  (X  -  I) 


(66a) 


Generally  in  a  supersonic  flow  ga  is  very  large>  so  that 
decreases  very  rapid]^  with  V,  therefore  expression  (66a)  can  he  written 
in  the  following  approximate  form: 


So  =  logp  ga 


\l2  -  1 


The  Integral  of  the  right-hand  member  has  been  calculated  in  the  case  with 
fomiula  ( h2. )  and  lias  tlie  valua 


d|  =  -p  ~ 

Vn  ^00 


where  x  is  now  the  distance  in  absolute  value  from  the  leading  edge 
Hence  formula  (66)  becomes 


=  - 


d\fn  a 


(1  +  kM<„2)  _  (ga)'-^  \l^  logg  Ig 


{2-1 
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The  perturbation  increases  In  absolute  value  from  the  leading  edge 

to  the  trailing  edge  with  an  average  gradient 


dx 


1  +  kM  2) 
00  • 


loge 


\/2  -  1 

2 


4ga  U 


(68) 


It  is  seen  that  in  a  supersonic  flow  the  heat  sources  will  produce  a 
decrease  in  the  velocity  of  the  flow  which  passes  over  them. 


NONLINEAE  COOLING  OF  A  HEATED  BODY  AT  SUEERSONIC  SPEED 
Nonlinear  Cooling  of  a  Flat  Plate 


The  role  of  the  heat  sources  along  a  body  in  a  compressible  flow  can 
be  imagined  to  heat  a  thin  layer  of  air  which  pertiirbs  the  motion  of  the 
main  flow,  just  as  if  the  body  were  deformed  instead  of  being  heated. 

The  average  velocity  which  passes  along  such  a  deformed  body  is  then* an 
"effective  velocity"  defined  by 


or 


U 


eff 


U 


+ 


dx 


Ueff/U  =  1 


+ 


1 

U 


d$^ 

dx 


where  dO^/dx  is  given  by  formula  (55)  in  a  subsonic  flow,  and  by 

formula  (68)  in  a  supersonic  flow.  For  the  sake  of  convenience  of 
writing,  introduce 


g  ^  ^eff  ~  ^  ^  1  ^ 

U  U  dx 

as  the  correction  coefficient  of  the  effective  velocity 


(69) 


Ueff/U  =  1  +  e 
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Introducing  this  effective  velocity  in  the  expression  of  heat 
delivery  (43)^ 


\j8nga{  1 


+  e/2)  +  1 


(70) 


where 


2-k 

"k-1 

g  =  1^00  hU 


In  problems  of  heat  transfer,  the  following  parameters  play  an 
important  part: 

(l)  The  Nusselt  number  N  as  introduced  previously  in  formula  (4o) 
is  defined  by 

Q  =  NXA  d^/d 

where  Q  is  the  amount  of  heat  delivered  in  \mit  time  by  the  unit  length 
of  an  immersed  body  across  an  area  A,  is  the  difference  between  the 

temperatvire  of  the  body  and  the.  isentropic  stagnation  temperature,  and 
d  is  a  representative  length.  For  a  flat  plate  of  length  4a,  N  is 


N  = 


Q 


(2)  The  Reynolds  number  R  =  pQUd/ri,  where  U  is  the  free-stream 
velocity,  Pq  is  the  isentropic  stagnation  density,  and  t]  is  the 
coefficient  of  viscosity 

(3)  The  Prandtl  number  Pr  =  2hT]/po,  where  2h  =  PqC-^/X 

(4)  The  Peclet  number  P  =  R  X  Pr  =  8hUa 

It  is  to  be  noted  that  if  the  Prandtl  number  is  a  constant,  the 
Peclet  number  is  proportional  to  the  Reynolds  number.  If  the  results 
are  to  be  expressed  in  terms  of  P,  the  quantity  8ga  which  occurs 
frequently  in  this  paper  becomes 

2-k 
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Hence  an  alternative  form  of  formula  (70)  for  the  heat  delivery  is 
obtained  as  follows: 


where  e  is  given  by  formulas  (55)  and  (68)  for  the  subsonic  and  the 
supersonic  speeds,  respectively,  as  follows; 

(a)  In  a  subsonic  speed: 


MI 

16  ''  1  - 


(b)  In  a  supersonic  speed: 

e  =  -i  (1  +  kM„2)(M^2  _  1)2  ^  ^ 

(72) 

By  substituting  formulas  (7I)  and  (72),  formula  (70a)  can  be  written 
also  in  the  following  form: 


2N  -  Cq^/ttP  +  -  (c2  logg  ^  “  C 


(70b) 


where  N  is  the  Nusselt  niimber,  P  is  the  Peclet  number  (it  is  equal 
to  the  Reynolds  number,  when  the  Prandtl  number  is  unity),  d  is  the 

difference  between  the  temperature  of  the  body  and  the  isentropic  stag¬ 
nation  temperature  Tq,  and  finally  Cg,  Ci,  Cg,  and  03  are  the 

following  functions  of  the  free-stream  Mach  number  M^  (see  figs.  3 
and  4) . 
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(2)  N  decreases  as  increases 

(3)  N  varies  linearly  with  \fF 

(4)  For  very  moderate  even  for  «  0,  formula  (70b) 

reduces  to 


2N  =  ^  +  1  +  ^ 


It  reduces  to  the  formula  of  King  (reference  5)  only  if  \/Tq  vanishes, 
(b)  Supersonic  flow: 

(1)  N  decreases  as  ■d^  increases.  The  rate  of  change  depends 

on  loge  /P,,  and  therefore  is  larger  in  absolute  value  than  the  rate  of 
change  in  a  subsonic  flowj  Q  does  not  vary  linearly  with  (fig*  5) 

(2)  N  decreases  as  Moo  Increases.  This  is  true  for  a  flat 

plate  with  a  sharp  nose.  It  will  be  seen  in  the  following  section  that 
for  a  body  with  a  blxmt  nose  different  aspects  will  appear 

(3)  N  does  not  vary  linearly  with  \fP)  because  the  temperature 

term  in  formula  (70b)  contains  a  function  logg  \fP.  However,,  the  devi¬ 

ation  from  the  linearity  is  small  owing  to  the  small  magnitude  of  -S^/Tq. 


Some  Remarks  on  Heat  Delivery  by  a  Body  With 
a  Blunt  Nose  in  Supersonic  Flow 

In  the  preceding  sections^  a  certain  continuous  distribution  of 
heat  sources  immersed  in  a  flow  of  velocity  U  has  been  considered.  By 
the  use  of  a  generalized  potential  theory  developed  for  a  nonadiabatic 
and  rotational  fluid,  the  action  of  the  heat  sources  on  the  flow  has  been 
studied,  and  the  heat  delivery  has  been  calciilated. 

In  order  to  show  in  the  theory  the  fundamental  action  of  heat  on 
the  flow,  and  to  obtain  the  characteristic  physical  features,  it  is 
believed  that  the  extraneous  action  of  the  shape  of  the  medium,  that  is, 
the  metallic  body  introducing  heat  to  the  fluid,  must  be  eliminated  in  ^ 
a  first  study. ^  Therefore,  in  the  present  calculations,  a  fine  plate  has 
been  chosen  which  does  not  offer  any  shape  perturbation  to  the  flow. 

This  circumstance  approaches  the  theoretical  distribution  of  heat  sources 
without  Interference  from  an  extraneous  body.  Fundamental  investigations 
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of  the  physical  phenomena  under  such  a  simplified  circumstance  may  serve 
as  a  basis  for  the  study  of  the  phenomena  under  other  more  complicated 
circumstances^  for  example^  the  action  of  heat  combined  with  the  action 
of  an  extraneous  body.  Hot-wire  anemometry  in  an  incompressible  fluid 
takes  usually  a  cylindrical  wire  for  the  introduction  of  heat  to  the 
fluid.  No  difficulty  has  ever  arisen  since  its  body  perturbation  is 
known  by  elementary  theory.  In  a  supersonic  flow,  the  body  perturbation 
is  more  difficult  to  calculate.  However,  for  the  study  of  the  heat 
delivery  an  approximate  estimation  can  be  made  by  starting  from  the 
following  considerations. 

A  blunt  body  like  a  cylinder  will  produce  a  detached  shock  wave 
which  is  strong  on  the  front  region  and  weak  on  the  two  branches .  Behind 
the  shock  wave,  the  flow  about  the  cylinder  is  subsonic  on  the  front 
region  and  supersonic  on  the  rest  of  the  region.  From  the  present  cal¬ 
culations  (see  section  Flux  Distribution  and  Heat  Delivery)  it  is  known 
that  the  air  in  the  front  region  gives  the  largest  contribution  in  the 
heat  delivery.  That  part  of  the  air,  heated  under  a  subsonic  state, 
passes  the  sonic  zone  and  leaves  the  cylinder  with  a  supersonic  velocity. 
While  the  detailed  picture  of  the  phenomena  has  to  be  studied  more  rig¬ 
orously,  a  rough  estimation  of  the  heat  delivery  can  already  be  obtained 
by  giving  to  that  contributing  part  of  the  air  an  effective  velocity 
after  a  Mach  wave,  and  an  effective  subsonic  Mach  number  Mp  which  can 
be  derived  from  the  unperturbed  Mach  number  Mo,  before  that  part  of 
the  shock  wave  which  is  approximately  normal.  The'  known  relation  between 
such  Mach  numbers  is  (see  fig.  4) 


1  + 


ilh) 


The  temperature  distribution  about  a  cylinder  is  given  by 


6^'^  4-  ^  ^  (75) 


This  is  the  general  form  of  equation  (23).  By  choosing  a  system  of 
coordinates  .  cp  and  t,  this  distribution  reduces  to  the  temperature 
distribution  in  a  uniform  stream  flowing  along  a  flat  plate,  as  given 
previously  by  equation  (23).  The  same  method  used  in  the  sections 
Temperat-ure  Distribution  and  Flux  Distribution  and  Heat  Delivery  can  be 
applied  here  in  order  to  calculate  the  temperature  distribution  and 
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finally  the  heat  delivery  for  a  cylinder.  Omitting  the  details,  the 
result  is 

Q  =  q/U  +  l)  (76) 

The  limits  of  variation  of  the  potential  function  are  0  and  3^.  In 
the  case  of  an  elliptic  cylinder  of  semiaxes  (a,b), 

3o  =  2(a  +  h)U 


As  pQ  -  4Ua  for  both  cases  of  a  circular  cylinder  of  radius  a  and  of 

a  plate  of  length  4a,  it  is  seen  that  the  expressions  of  the  heat 
delivery  for  both  cases  without  heat  perturbation  must  be  identical.  It 
is  also  to  be  noted  that,  for  =  0,  the  expression  of  the  heat  deliv¬ 
ery  reduces  to  that  found  by  King.  Further,  for  a  circular  cylinder. 


N  =  — - — 

=  4hUa 


(Pc)e„  = 


Hence  formula  (76)  becomes 


\/2 


(77) 


The  mass  of  air  on  the  front  part  of  the  cylinder^  heated  at  a 
subsonic  state  characterized  by  the  subsonic  Mach  number  will  leave 

the  cylinder  with  a  supersonic  velocity  which  is  approximately  the 
effective  velocity  after  a  Mach  wave 


Ueff  =  U(1  +  e) 
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where  e  is  given  by  expression  (72).  By  substituting  expression  (72) 
into  formula  ( 77) ^ 


^  4-  i  r-,  _  Inc.  »/p"  - 


Nc  =  V  Co\/^  +  7  Cl  -  ■^/Cg  loge  ^  -  03]^ 


where  the  C ' s  are  the  following  functions  of  the  Mach  numbers 


Co(Mi)  =  Pi 


Ci(Mi)  =  Pi"-^ 

Co(Mi) 


C3(“l)  = 


2  -  k 


loge  Poo  -  loge 


\/F-  1 
4 


Formula  (78)  remains  valid  for  a  subsonic  flow,  provided  the  values 
of  the  C-functions  are  taken  from  expressions  (73). 

The  amount  of  heat  Q  does  not  vary  with  (see  fig.  7)^  although 

the  deviation  from  linearity  is  not  so  important  as  for  the  case  of  a 
flat  plate. 

For  small  values  of  the  body  perturbation  becomes  predominant 

compared  vith  the  heat  perturbation.  Therefore  increases  with  Moo 

for  Mqo  >  1,  a  behavior  directly  opposite  to  that  for  a  flat  plate 
which  does  not  offer  any  body  perturbation  (see  fig.  8).  On  the  other 
hand^  for  large  values  of  ^Sv/^o  (see  fig.  9)^  the  body  perturbation 
becomes  less  important  compared  with  the  heat  perturbation.  Therefore 
Nc  decreases  with  Moo  for  Moo  >  1^  a  behavior  comparable  with  that  for 

a  flat  plate  (cf.  fig.  6). 


From  the  above  formula  for  the  heat  delivery,  it  is  seen  that  Q 
is  proportional  to  the  coefficient  or  x  Re.  In  the  present 

theoretical  investigation,  this  coefficient  has  been  assumed  as  independ¬ 
ent  of  the  temperature  gradient  in  the  neighborhood  of  the  wire,  and 
the  viscosity  has  been  assumed  negligible.  In  view  of  this  simplification 


NACA  TN  2436 


49 


and  of  the  uncertainty  attached  to  the  value  of  X,  the  theoretical 
value  of  the  coefficient,  calculated  for  a  constant  temperature  of  170  C, 
seems  too  high  compared  with  the  experimental  ones.  By  comparing  the 
formula  of  King  with  the  experiments  on  cylindrical  hot-wires  at  zero 
Mach  number,  it  is  seen  that  the  ratio  of  experimental  coefficient  to 
theoretical  coefficient  is  in  the  range  of  1,5  to  1.7.  These  consider¬ 
ations  have  to  he  taken  into  account  in  the  numerical  interpretations  of 
the  results. 


Some  Concluding  Remarks  on  Heat  Flow  and 
Some  Notes  on  Effects  of  Viscosity 

In  this  section  various  details  of  the  reasoning  are  reviewed,  and 
the  effects  of  viscosity  are  studied.  This  will  fvirnish  at  the  same  time 
an  opportunity  of  pointing  out  briefly  the  various  principles  and  hypoth¬ 
eses  underlying  the  theory,  and  the  difficulties  which  herewith  are 
connected. 

In  the  present  paper,  it  was  Intended  to  elucidate  some  fundamental 
theoretical  features  of  a  heat  flow,  namely, 

(1)  The  generalized  potential  theory  for  a  flow  with  rotation  and 

variable  heat  energy 

(2)  The  distribution  of  heat  energy  fvirnished  by  a  heat  source 

(3)  The  distribution  of  heat  fltix 

(4)  The  amount  of  heat  delivery 

(5)  The  perturbations  of  the  flow  by  heat 

A  simplified  model  was  chosen  for  the  study;  it  consisted  of  a 
straight  array  of  continuous  heat  soiorces,  placed  parallel  to  a  stream 
of  uniform  velocity.  The  heat  sources  on  the  array  were  of  such  an 
intensity  distribution  as  to  furnish  a  prescribed  surface  temperature, 
for  instance,  a  uniform  temperature.  In  order  to  fix  the  idea,  such  an 
array  was  called  an  "infinitely  thin  flat  plate,"  which  has  the  only 
significance  of  indicating  the  location  of  the  heat  sources  giving  a 
constant  surface  temperature,  without  offering  any  body  perturbation  to 
the  flow,  for  example,  in  the  form  of  a  bovindary  layer. 

The  theory  of  generalized  potential  flow  has  the  advantage  of  per¬ 
mitting  the  use  of  the  potential-streamline  coordinates,  and  hence  the 
application  to  more  complicated  shapes  of  heat  source  distribution,  other 
than  the  above  straight  array.  In  fact  the  use  of  such  coordinates  can 
transform  a  complicated  distribution  into  the  straight  array,  which  is 
therefore  the  most  fundamental  distribution. 
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The  simplified  model  has  served  to  illustrate  certain  physical 
features.  However,  it  may,  of  course,  deviate  appreciably  from  an  actual 
body  supplying  heat  to  the  flow.  For  example,  the  flow  in  the  boundary 
layer  of  the  body  should  be  actually  retarded,  and  therefore  the  heat 
delivery  should  be  less  intense.  But  before  going  to  such  problems  of 
a  nonuniform  stream,  it  was  necessary  first  to  study  the  more  fundamental 
problem  of  a  uniform  stream,  which  gives  a  clearer  analysis  of  the  pre¬ 
dominant  physical  characteristics.  The  heat  delivery  by  a  blunt  body, 
for  example,  a  cylinder,  gives  rise  to  additional  difficulties  which  were 
not  investigated,  because  the  nonuniform  flow  behind  a  curved  detached 
shock  wave,  as  produced  by  the  blunt  body,  is  in  itself  also  a  difficult 
problem,  even  without  heat  introduction,  and  must  be  left  out  of  the 
scope  of  the  investigation  at  the  present  time.  In  that  respect  only 
some  phenomena  behind  such  a  wave  have  been  roughly  stipulated  in  the 
preceding  section,  illustrating  the  difficulties  which  arise,  rather  than 
their  explanations. 

The  theory  of  the  heat  flow  ( nonadiabatic,  rotational,  and  with 
variable  total  energy),  in  connection  with  the  above  enumerated  aspects, 
was  elucidated  in  the  preceding  sections  by  neglecting  the  effects  of 
viscosity,  a  procedure  generally  adopted  in  the  theories  of  such  a  com¬ 
plicated  flow.  In  the  following  lines  the  effects  of  viscosity  will  be 
estimated. 

For  a  viscous  flow,  the  equation  of  energy  is  (see  reference  6) 


DT 

Dt 


+ 


where 


Dt 


—  +  w  -  grad 
ot 


$  is  the  dissipation  function  is  used  for  potential  perturbation, 
and  the  symbol  ^  is  therefore  chosen  for  the  dissipation  function) 
and  all  other  symbols  have  been  explained  in  the  first  section  of  the 
present  paper. 

By  substituting  for  p  from  the  Navier-Stokes  equation 


P  ^  =  -grad  p  + 
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where  t]  is  the  coefficient  of  viscosity. 


*  f )  -  I?  =  "  *) 


For  a  stationary  flow^ 


pw  •  grad  c^T  + 

\  P  2 


XV^T  +  ^TiwV^w  +  (79) 


This  partial  differential  equation  determines  the  spatial  distribution 
of  the  temperature  T^  or_,  following  equation  (l)^  the  spatial  distrib¬ 
ution  of  '3(x^y)^  that  is,  the  additional  temperature  furnished  by  the 
heat  source.  As  was  done  in  the  section  Small  Perturbations,  small 
perturbations  will  be  considered  by  keeping  only  terms  in  •6  of  the 
first  order,  and  hence  it  is  possible  to  transform  equation  (79)-  into  a 
linearized  partial  differential  equation  for  <9.  From  the  equation  of 
continuity  (5)^  and  from  the  relation  of  the  density-temperature  ratio, 
it  is  observed  that  Sw^/^Xj  is  of  the  same  order  as  ^ hence 

the  second-order  terms 


Sw.  Swi. 
o  X  j  u  ^  2 


=  1. 


2,  3) 


are  negligible. 

The  difference  between  equation  (79)  and  the  energy  equation  (6) 
used  formerly  lies  in  the  presence  of  the  last  expression  within  paren¬ 
theses,  which  represents  the  effects  of  viscosity. 

Now, 
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The  last  term  of  the  right-hand  member  is  a  second-order  term  in  t3  and 
is  therefore  negligible.  So  is  $  too.  Hence  by  neglecting  second- 
order  terms  in  the  last  expression  between  parentheses  of  equation  (79)^ 


pw  •  grad  fcpT  +  +  T]V^  ~ 


=  ^  +  Ip 

^  X  2  ; 


If  the  Prandtl  numher  is  written  as  Pr  = 


V^fc^T  +  — ^  ^  V  .  grad  fcpT  + 


The  density  p  can  be  transformed  into  T^  and  the  temperature  T  can 
be  further  transformed  into  using  the  formulas  of  transformations  (2) 
and  (3)*  By  leaving  out  second-order  terms,  after  transformation  there 
is  obtained 


= 


2g' 


where 


(80) 


_ 

f  ^"1  t.TT  -1 

g  =  hco 


and  g,  h,  and  U  are  defined  as  usual  in  the  first  section  of 

the  paper.  Equation  (8o)  is  the  partial  differential  equation  for  the 
spatia-1  distribution  of  by  taking  into  account  the  effects  of  vis¬ 
cosity.  It  is  seen  that  equation  (8o)  is  of  the  same  form  as  the  non- 
viscous  distribution  represented  by  equation  (23)^  except  the  coeffi¬ 
cient  2g  in  equation  (23)  is  replaced  by  2g'  in  equation  (8o). 

The  spatial  distribution  of  the  temperature  '<3,  the  flux  distrib¬ 
ution,  and  the  amount  of  heat  delivery  follow  from  equation  (80),  exactly 
in  the  same  forms  as  formulas  (29)^  (STa),  and  (43)  with  the  new  coeffi¬ 
cient  2g'  as  follows: 
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-d  =  — e®  (ff'r)u 
2nX  o''8 


q(x)  =  q(0)e^®'^(l  -  erf  \/2g^) 


with  q(0)  =  2X\  2g'/n^  (  according  to  formula  (4l)) 


N 


=  ^(fSrrg'a  +  l) 


(81) 


It  is  seen  that  the  viscous  distributions  of  i3  and  q  are  "rounded 
down^"  that  is^  not  so  sharp  as  the  nonviscous  distributions,  and  that 
the  heat  delivery  is  diminished  by  the  viscosity.  The  effects  of  the 
viscosity,  which  increases  with  the  free-stream  Mach  number,  are  thus 
to  modify  the  magnitude  of  the  physical  quantities  considered  above, 
without  changing  the  forms  of  their  differential  equations  and  their 
formulas . 

From  formula  (8l)  for  the  heat  delivery,  it  is  seen  that  N  is 
proportional  to  and  hence  to  Now  the  horizontal  velocity  U 

is  subject  to  perturbations  by  the  heat  transfer.  Those  pertiirbations, 
of  the  order  of  ■fl,  were  assumed  small,  and  as  a  first  approximation  they 
are  neglected  in  the  elucidation  of  formulas  (43)  and  (8l)  for  the  heat 
delivery  in  a  nonviscous  flow  and  a  viscous  flow,  respectively.  Later, 
as  a  second  approximation,  the  flow  pertiirbations  by  heat  were  calculated 
in  the  sections  explaining  flow  perturbations  by  heat  sources  at  subsonic 
and  supersonic  speed,  by  neglecting  the  viscosity.  It  is  an  easy  matter 
to  substitute  the  perturbed  velocity  in  the  formula  of  unperturbed  heat 
delivery  in  order  to  find  the  perturbed  heat  delivery.  Probably  the 
viscosity  will  again  have  certain  effects  on  the  flow  perturbations. 

But,  as  the  flow  perturbations  play  only  a  second-order  part  in  the  heat 
delivery,  those  effects  will  not  be  studied  further  here. 

Curves  are  plotted  in  figures  for  the  case  of  a  nonviscous  flow, 
illustrating  the  general  properties  of  a  heat  flow.  The  corresponding 
curves  for  a  viscous  flow  are  not  added,  in  order  not  to  make  the  figures 
too  cumbersome. 

The  effects  of  viscosity  will  reduce  the  slope  of  the  curves  for 
heat  delivery  in  figures  6,  8,  and  9,  without  changing  the  intercepts 
with  the  vertical  axis.  The  slope  of  the  curves  in  figure  6  will  have 
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to  be  multiplied  by  a  factor  is,  where  Cj  is 

given  by  figures  3  and  4.  The  slope  of  the  curves  in  figures  8  and  9 
will  have  to  be  multiplied  by  the  factor 


-1/2  ^ 


where  is  given  by  figure  4  in  terms  of  Moo* 
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Figure  4.-  Mach  number  functions  at  a  supersonic  speed. 
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Figure  5- -  Nonlinearity  of  heat  delivery  at  a  subsonic  and  a  supersonic 
speed  in  the  case  of  a  flat  plate.  Broken  lines  represent  corre¬ 
sponding  linear  relations  obtained  by  neglecting  heat  perturbations. 
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linear  relations  obtained  by  neglecting  heat  perturbation.  Peclet 

number  for  a  cylinder. 
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Figure  8.-  Eelation  iDetween  heat  deliveiy^  •  the  Peclet  number,  Mach 
number,  and  temperature  difference  in  the  case  of  a  cylinder  for 
=  0.05.  The  relation  for  Moo  =  1  and  -^v/^o  ^ 

plotted,  as  it  is  not  covered  by  the  elliptic  and  hyperbolic  equations 
of  perturbations. 
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The  relation  for  Moo  =  1  and  t3v/To  f  0  is  not  plotted,  as  it  is 
not  covered  hy  the  elliptic  and  hyperbolic  equations  of  perturbations. 
For  large  Mach  numbers  the  theory  does  not  predict  variations  of 
Nusselt  numbers  at  low  Peclet  numbers. 
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